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ABSTRACT 

We present the initial results of a 40 night contiguous ground-based campaign 
of time series photometric observations of a 1.39 deg^ field located within the 
NASA Kepler mission field of view. The goal of this pre-launch survey was 
to search for transiting extrasolar planets and to provide independent variability 
information of stellar sources. We have gathered a data set containing light curves 
of 54,687 stars from which we have created a statistical sub-sample of 13,786 stars 
between 14 < r < 18.5 and have statistically examined each light curve to test 
for variability. We present a summary of our preliminary photometric findings 
including the overall level and content of stellar variability in this portion of the 
Kepler field and give some examples of unusual variable stars found within. We 
present a preliminary catalog of 2,457 candidate variable stars, of which 776 show 
signs of periodicity. We also present three potential exoplanet candidates, all of 
which should be observable in detail by the Kepler mission. 
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Introduction 



Stellar variability studies provide critical access to a number of astronomically significant 
properties including rotation rates, eclipsing binaries, and pulsations. Aside from offering 
insight into the nature of the stars themselves, statistical studies conducted on large samples 
of field and cluster stars broaden our understanding of stellar evolution and can assist in the 
search for exoplanets. Variability surveys that endeavor to find variations with intermediate 
periods of days to weeks studying objects such as short period eclipsing binaries and pl anetary 
trans i ts, benefit from dedica ted long-term, high cadence, dedicated observations (IHowell 
20081 : Ivon Braun et al.l l2009bl and references therein) . 



NASA's Kepler mission ( iBorucki et al.ll2010af ). which was launched in March 2009, is 
conducting a transit search in Cygnus with the goal of finding Earth-like planets orbiting in 
the habitable zones of sun-like stars. The signatures of transits of this nature would have 
relatively small depths, making the inherent variability of the host star even more relevant. 
The success of this mission partially depends on an accurate characterization of the stellar 
variability in the field. To that end, the Burrell-Optical-Kepler-Survey (BOKS) was designed 
to determine the level and type of stellar variability in a small (~ 1%) subsection of the Ke- 
pler field. As an added goal, we can assess the frequency of close-in Jovian-type planets 
(the so-called "Hot Jupiters" ) in the same field, and allow for a comparison of ground-based 
and the Kepler based transit surveys. A number of other Hot Jup iters have been discov- 



ered in the Kepler field prior to l aunch from ground-based surveys ( jO'Donovan et al.ll2006 



Pal et all 120081: iBakos et all 12010 ) and Kepier itse l 



exoplanets ([ Borucki et al. 



2010b 



Jenkins et al.w2010i : lBorucki et al. 



Koch et al. 



2010c 



2010 



has already discovered many additional 



Steffen et al. 



Dunham et al. 2010; Latham et al. 2010 



20101 ) . Further comparisons of ground- 



based and space-based transit candidates would be extremely beneficial due to the high qual- 
ity lightcurves that Kepler can provide. In particular, since the Kepl er mission must mak e 
numerous selection cuts in order to achieve their mission objectives ( iBatalha et al.l 120101 ). 
there is a tendency to avoid fainter target stars that may have detectable Hot Jupiters, 
but not be suitable for Earth-sized transit searches. By identifying additional Hot Jupiter 
candidates, and then having Kepler undertake follow-up observations, we may be able to 
study detailed properties of these systems and characterize them in exquisite detail. This 
has already been done for one pre-launch Hot Jupiter candidate, HAT-P-7, and there are 
indic ations that the ex trasolar planet in this system is gravitationally distorting the host 
star dWelsh et al.lboiol ). 



This paper introduces the properties of BOKS, and gives a summary of the data reduc- 
tion and analysis of the survey. In [|2]we present our observing strategy and a summary of 
our observations. We outline our data reduction techniques and discuss the observational 
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window function of our survey in ^ We present the object detection, photometry and as- 
trometry in ^ Finally, in ||5] we discuss the initial results of our variability survey and of 
the search for exoplanets in the BOKS field, specifically those of the "Hot Jupiter" variety. 
We review our conclusions in |JHi There are a number of other scientific projects planned 
for the BOKS data, such as the comparison of stellar field to cluster variability, cataloging 
the many variable stars found in this survey, and searching for any moving objects. These 
results will be discussed in future papers. We plan to submit all of the BOK S data to the 
NASA/IPAC/NExScI Star and Exoplanet Databas^ ( von Braun et al. 2009a ). where it can 
be of service to the entire astronomical community 



2. Observations 



In order to maximize the scientific benefits from this survey, we chose our scientific 
field under a number of constraints. First, besides being located within the Kepler field our 
field ideally should have a large number of stars, as this will improve our chances to find 
extrasolar planets and find other interesting objects. Second, given that our ground-based 
imager has relatively large pixels in angular size, the field must not be so close to the Galactic 
plane that photometric crowding would be a major factor. Third, in or der to compare ou r 
data against num erous stellar cluster var iability surv eys, such as UStAPS(|Hood et al.ll2005[) . 
EXPLORE/Q C jvon Braun et all 120051 ). PISCES JMochejska et"ld1 boOfil ). and STEPSS 
( Burke et al. 20061 ). we decided on a field that had both field stars and an open cluster 
within it, so that we could compare the variability properties of both stellar populations 
simultaneously. In order to find the optimal field, we first pre-imaged a number of candi- 
date fields on 2006 April 24. After visual inspection of all of the pre-images, we made a 
determination of the field that best matched our conflicting criteria. 

Our final selected target field of view covered 1.39 deg^ in the constellation of Cygnus and 
was centered on the open cluster NGC 6811 (RA= 19^37'"17^ Dec=+46°'23'^18"; WEBDA@). 
Our field is completely contained within the Kepler field, and the BOKS field is located on 
channels 63, 47, 23, and 39 of the Kepler imager in the Spring, Summer, Fall, and Winter 
seasons respectively §. Our photometric survey began on 1 September 2006 and ended on 
10 October 2006, consisting of 40 nights in all. We observed with the Case Western Reserve 



^located at: |http: / /nsted.ipac.caltech.eduT] 



^The WEBDA database, developed by J.-C. Mermilliod, can be found at http://www.univie.ac.at/webda/ 
^Full Frame Images (FFIs) of the Kepler field for each observing season can be found at 



http: / / archive.stsci.edu/kepler/fH_display.php 
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University 0.61m Burrell Schmidt telescope (hereafter the Burrell), located at Kitt Peak 
National Observatory. One advantage of a dedicated observatory for this sur vey is that 



we could observe for a large number of consecutive nights. Many researchers (jPont et al. 



20061 : iBeatty fc Gaudil 120081 : Ivon Braun et al.ll2009bl ) have shown that the total duration of 
a photometric survey is crucial in maximizing transit detection efficiency in the presence of 
statistically correlated ("red") noise. 

The imager used for this survey was a SITe back- illuminate d 2kx4k CCD with 15 micron 
(1.45 arcsecond) pixels, run by a Leach version 2 controller (ILeach et al.lll998l ). and two 
output amplifiers. The long axis of the CCD was oriented East/West. We observed primarily 
in the SDSS r-band filter, but we also obtained occasional Johnson l^-band images of the field 
of view in order to obtain two-color information and to allow for cross-comparison between 
our data and other photometric catalogs. In particular, we compared our photometry to 
that found from the Kepler Input Catalog (Kicfl 

Observations were ongoing during any weather conditions where stars were visible on 
the sky, and it was safe to operate the telescope. As a result, the BOKS data have large 
variations in seeing, transparency, and night sky brightness. Of the 40 nights of observing, 
thirteen were completely lost to weather, leaving 27 nights of potential data. Nearly all of 
our r-band and V integrations were 180 seconds in duration, with the exception of 31 images 
taken on night 14 that had exposure times of 300 seconds. The CCD readout time was 45 
seconds in length. A total of 1,924 r images and 10 V images were taken over the entire run. 
The SITe CCD gain was fixed at 2 electrons per ADU, each pixel had a full well capacity of 
at least 100,000 electrons, and the read noise was 12 electrons. We note here that while the 
CCD pixel scale is fairly large compared to typical CCD imagers, studies (Feldmeier et al. 
2011, in preparation) have shown that milli-magnitude relative photometry is possible, even 
with such large pixels. Table [T] provides a summary log of our observations and Figure [T] 
shows a graphical representation of the number of exposures throughout the survey. An 
r-band exposure of BOKS, created from co-adding the first 24 images from our survey is 
shown in Figure [2l 

With the help of the American Association of Variable Star Observers (AAVSO), we 
also arranged to have bright variable stars photometrically monitored in the field at the same 
time as BOKS was un derway. The prelirn inary results from this independent photometric 
survey are discussed in iHenden et al.l (120061 ) , and a more careful comparison will be discussed 



"^Version 10 of the KIC is available at: http://archive.stsci.edu/kepler/kiclO/search.php In some cases, 
we used the 7th and 8th version of the KIC for steps in our analysis. These will be referred to as KIC78 in 
the text. 
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in a future papei§|. 



3. Data Reduction 

Although the data reduction of BOKS is relatively straightforward, the large number 
of images and the need to ensure highly precise relative photometry demands some careful 
attention. We therefore began our CCD reductions as follows. Since our imaging observations 
were obtained using the dual amplifier mode with the SITe CCD, we first combined the two 
amplifier readouts into single images using IRAF'fl mscred.mkmsc task. The resulting images 
were then merged, trimmed, and overscan subtracted using the ccdred.ccdproc task. 

We obtained approximately ten bias and five twilight flatfield frames per night, which 
were used on the respective night's images using IRAF's ccdred package after they were 
checked for unwanted features such as bright stars in the flatfields or amplifier noise in the 
bias frames. If any such features were present, the corresponding bias or flatfield frames were 
discarded. We created nightly master bias frames, but due to the presence of dust grains on 
the dewar window and/or filter, that changed positions between nights, we did not create a 
master fiat for the entire run. 

After the data had been processed, we inspected the entire data set for quality issues. 
We determined basic parameters of each image such as the median sky level and the median 
seeing in order to verify that our data was suitable for use. To determine the median sky 
level, we used the IRAF task imstatistics in iterative mode. The median seeing of each image 
was calculated by applying the imexam task on 150 bright, unsaturated stars on each image. 
Each star was fit using a Gaussian function, and the median of all of the derived FWHM 
values was adopted for the median seeing for each image. 

After applying the photometric zero point (§4.3), we present the median sky brightness 
and median seeing for BOKS in Figure El As would be expected from a telescope run of this 
length, these properties varied substantially as differing lunar phases and weather patterns 
occurred. After this analysis, we removed night 37's data due to very high sky levels by 
the nearly full moon. Next, we visually inspected each image for quality issues. From this 



^The AAVSO NGC 6811 campaign information can be found at 
http: / / www.aavso.org/news/ngc681 l.shtm] 

^IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc, under cooperative agreement with the National Science 
Foundation. 



- 6 - 



process, we found that night 12 suffered from condensation on the CCD dewar window, 
leaving a circular distortion feature on approximately 17% of the area of each image. Night 
19 suffered from CCD electronics issues on one of the two CCD amplifiers. In both of these 
cases, many of the stars remain unaffected on each image. However, to be conservative in 
this initial study, we have removed the entire night's data from further consideration. After 
our initial light curve analysis (discussed in |J5]), we found that two additional nights, nights 
13 and 23, had extremely variable clouds (with such a large angular field, even the ensemble 
photometry method discussed in §l]can fail if the clouds are variable enough; Feldmeier et al. 
2011, in preparation). Although some of the exposures on these nights should be acceptable 
for our scientific goals, we again chose to be conservative and removed the entire night's data 
from consideration. This left us with imaging data from 22 different nights over the span 
of the survey and a total of 1,565 r-band images that could be used for potential variability 
studies. 

For a planetary transit survey, understanding the observational window function gives 
crucial insight into the survey completeness and sensitivity to periodic variability. Normally, 
the window function is defined as the probability that a planetary transit is detected in a 



given data set, as a function of planetary orbital period (jvon Braun et al.ll2009bl ). For BOKS, 
we calculated the approximate observable window function in the following manner: we 
simulated planetary transits from periods ranging from 0.5 to 30 days and divided each period 
into 10,000 phases. Given the starting and ending times of the 22 nights of observations, 
we then determined the probability of observing a transit with that period over all phases. 
The results of this analysis are plotted in Figure Si As with all time-limited photometric 
surveys, we are most sensitive to short period transits, and our ability to detect transits 
decreases with transit period. Since our survey is ground-based, the characteristic aliasing 
period of integer days is also strongly present in our data. We should note that this window 
function deals with temporal sampling only and likely to be an overestimate: it does not 
take into account the effects of differing transparency, seeing, and sky brightness on the 
detectability of transits. It also does not take into account the effects of differing stellar and 
planetary radii on the detectability of transits. Finally, the effects of statistically correlated 
"red noise," which are likely to be significant, are not included in this calculation. We plan 
to perform extensive Monte Carlo simulations on these effects, which will be presented in 
a future paper. However, we note that BOKS has a significant advantage over many other 
ground-based transit surveys: if even one transit appears in our survey, we could, in principle, 
verify it with Kepler follow-up observations. 



-7- 



4. Object Detection, Photometry, and Astrometry 

With our good quality data set finalized, we next focussed our attention on finding all 
stellar sources in the BOKS field, and determining their magnitudes and positions throughout 
the survey period. 

4.1. Object Detection and Coordinate Transformations 

We next created a master image by combining six individual images from night 25 using 
the imcombine task. This master image was used to ensure that we detect all of the stellar 
sources in the frame and remove the possibility of radiation events contaminating our source 
catalog. We next created the master list of stars (point sources) by running the daofind task 
within IRAF on the master image. We chose a threshold value for point source detection 
of five times the standard deviation of the sky background. Given the large pixel scale of 
our data, we adjusted the sharphi value to 0.9 rather than the 1.2 that is normally adopted. 
All other data-independent parameters were left at their default values. We found a total 
of 56,354 sources that matched the daofind criteria we adopted. This initial list was then 
manually inspected to ensure that non-stellar objects were not included, such as diffraction 
spikes, radiation events, or objects on the extreme edges of the frame. This left a total of 
54,687 objects for further study. 

With the master list of coordinates determined, we then needed to re-identify each source 
on every frame of our survey. Rather than shift each image, which would lead to unacceptable 
uncertainties in the magnitudes due to interpolation, we instead re-determined the source 
coordinates for each individual image. To do this, we used a high quality image from the 
middle of night 25 as our positional reference image. To match stars in individual images to 
those found on the reference image, we split each image into eight rectangular subsections 
(512x512 pixels). For each subsection of each image, we summed a group of rows and group 
of columns to create a pair of one dimensional arrays that contained the peak counts at the 
row or column location of each star. We then used the Fourier transform between pairs of 
these arrays taken from different images to locate translational shifts between each subsection 
of the images. These shifts were applied to the master coordinate list to locate each star on 
individual images. The Burrell tracked fairly well overall but in periods of cloudy conditions, 
when the autoguider was unable to hold the tracking steady, the spatial shifts could be up 
to 1-2 arcminutes. By dividing the image into subsections, we were able to accommodate 
most of the magnification and rotational differences between images in our final astrometric 
coordinate solution. 
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4.2. Aperture Photometry and Ensemble Correction 

We performed aper ture photometry using the phot task within IRAF's noao.digiphot.daophot 



package (IStetsod 119871 ) . We experimented with five seperate photometry apertures, from 
very small radii (1 pixel) to large radii (5 pixels), to span the range of stellar brightness and 
seeing changes over the observed field of view. After inspection of the output light curves, 
we selected two distinct aperture values for our work (3 and 4 pixels) based on their small 
level of scatter around the median magnitude of a representative light curve sample. These 
two aperture values, which correspond to 4.35 arcseconds and 5.8 arcseconds, with a sky 
annulus of 11.6-33 arcseconds, gave good results and we used the smaller aperture for our 
final light curve set. 

Once the raw photometry files were created, we then needed to account for the effects 
of seeing, transparency, and airmassfor each star. To do this, we adopted a local ensemble 
photometry technique, where we used a local set of bright stars that are photometrically 
co nstant to determine chan ges in these param e ters. The algorithms are discussed in detail 



by lEverett et al.l (120021 ) and lEverett fc Howell! (120011 ): we will briefly outline them here. 



First, we divided each image up into 8x4 square regions with a size of 500 x 500 pixels 
(corresponding to 725 arcseconds on the side). This size was chosen to allow a sufficient 
number of ensemble stars to be present in the individual regions and at the same time, to 
optimize our sampling of positional dependence of photometric effects such as variable point 
spread functions, color terms, or focus gradients. For a star to be an acceptable ensemble 
star, it must fulflU certain quality criteria. Speciflcally: 1) the star must be present in all 
frames, 2) it must have a photometrically constant light curve (x^ < 3.), 3) it must be bright 
enough that photon noise is negligible (average flux must be greater than 50,000 ADU, 
corresponding to a S/N of 316 or better), and 4) have no close-by stellar companions that 
would interfere with the light curve in poor seeing conditions (no stars within 5 pixels that 
are within 5 magnitudes of the ensemble candidate's magnitude). We created an initial list 
of ensemble stars by conducting an automated search through the stars in each region. After 
this automated preselection, the ensemble stars were inspected by eye to eliminate stars that 
appeared to have signs of residual variability compared to the remainder of the ensemble 
stars. A total of 688 stars with r magnitudes varying between 14.3 and 16.6 were used in 
the flnal ensembles. 

After the flnal list of ensemble stars was created, the relative photometry procedure was 
rerun using only the cleaned sample of stars. If a region had fewer than ten ensemble stars, 
we combined it with a neighboring one to create a larger region. The exact calculation of 
relative photomet ric offsets for t he in dividual regions was performed by a custom written 



routine based on lEverett et al.l (120021 ). An example of these offsets for a single region is 
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plotted in Figure |5l Due to the effects of airmass, seeing, and transparency, these offset 
values vary significantly over the survey length. 

Due to the differing positions of each individual image, various stars have differing 
numbers of observations, with objects at the edges of the fields having fewer observations 
than object near the center. To ensure a high quality set of light curves for study, we focus 
exclusively on light curves that have at least 1,000 photometric measurements. This left 
32,806 sources for further analysis. 

An unfortunate issue we found during the data reduction was the fact that the start 
time (HJD) for any given exposure listed in the image header was not sufficiently accurate 
for our purposes. The internal clock on the data recording cornputer, which was a Microsoft 
Windows PC running the Voodoo image acquisition softwarqj at the telescope was found 
to be imprecise and could not be used alone for exact timing purposes as we found that it 
drifted by up to several minutes over the course of a single night. As a result, we used the file 
creation times (recorded by a different clock on an internet time controlled Linux machine) 
as HJD "start time" information in our image headers. From some experimentation, we 
found that the header start time recorded is approximately 93 seconds after the true mid- 
exposure time in most cases. From comparing consecutive exposures throughout the survey, 
the uncertainty in this correction is approximately one second. Consequently, the times 
we recorded in our light curves will be precise relative to each other, but our time zero 
point is not tied to UT or any other absolute time system to within several minutes. In 
the future, we plan to correct this by correlating our observations against the AAVSO and 
Kepler observations, which should allow us to reduce any time offsets. 



4.3. Photometric Zero Point 

The transformation from differential instrumental magnitudes to SDSS r magnitudes 
for our stars was performed by comparing the magnitudes of the ensemble stars in each 
subregion to the corresponding r stellar magnitudes listed in the Kepler Input Catalog (KIC). 
The scatter in this comparison was approximately 0.1 magnitudes, which we take as our 
systematic magnitude uncertainty for the BOKS survey. Given that we have minimal color 
information in our observations and our survey is primarily interested in searching for stellar 
variability, this amount of uncertainty is acceptable for our needs. For any individual object, 
however, one can derive a more accurate magnitude by directly comparing the object to the 



''located at ,http: / /www.astro-cam.com/ 
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KIC values, which have a mean r systematic uncertainty of less than 0.015 magnitude^. 



4.4. Astrometry 



We performed astrometry using the USNO-Bl.O catalog ( iMonet et al.ll2003[ l and wc- 
stool^. In order to eliminate the effects of field rotation and distortion, we performed 
astrometry on the individual regions as described above, rather than on the field as a whole. 
Our internal uncertainties on our coordinates were approximately ± 0.5"right ascension and 
± 0.3"declination. As an independent check on our astrometry, we co mpared our final astro - 
metric catalog to that found from the 2MASS Point Source Catalog (jSkrutskie et al.ll2006l ). 
Figure [6] shows the astrometric comparison between our BOKS astrometric solutions and 
those listed for common stars in the 2MASS catalog. The one-sigma coordinate offset is 
near 0.4"both RA and DEC with a small, but clear, asymmetry across the field. Given that 
the pixel size of the Burrell is 1.45"and that the wide field of the Burrell produces differential 
focus and refraction effects, this agreement is quite reasonable. 



5. Light Curve Analysis 

With the light curves established for each star, we then began the search for photometric 
variability. To characterize the possible variable nature of our light curve sample and to 
search for possible exoplanet transits, we utilized OPTICSTAT, a custom-written statistical 
analysis package created by M. Everett (discussed in Howell et al. 2005). OPTICSTAT 
returns several statistics related to stellar variability, including the reduced (x^)? the 
standard deviation, the probability that the star varies periodically, and the most likely 
period. Prior to our statistical light curve analysis, we removed the effects of obvious cosmic 
rays which will artificially increase the apparent flux as follows. If there were one or two 
consecutive points in the light curve that deviated from the mean magnitude by more than 
3.5 times the standard deviation, those points were rejected during statistical tests. This 
might remove some signs of ultra-short variability, but since the exposures are 180 seconds 
in length, and the readout time of the CCD is 45 seconds, the total scientific impact should 
be minimal. 

To test each light curve for general variability, we fitted the light curve with its mean 



^photometric uncertainties for the KIC can be found at: http: / /www. cfa.harvard.edu/kepler/kic/kicindex. html 



http : //tdc-www . harvard . edu/sof tware/wcstools/ 
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flux and then calculated the probability that the reduced xl statistic shows the data to 
be inconsistent with this mean flux. This test can easily fail, however, in the presence 
of systematic errors or uncertainties in the calculated magnitude errors we assign to the 
individual data points. To compensate, we adopted a more conservative threshold for the 
xl probability than our formal errors would dictate. When applied to the full light curve, 
we have adopted the criteria that point sources that have xl ^ ^ cire variable sources. 



5.1. Photometric Uncertainties 

In order to make a proper determination of variability for each star, we must deter- 
mine the random and systematic uncertainties for our target stars. By splitting the entire 
BOKS field of view into 32 smaller sub-sections we were able to remove the vast major- 
ity of systematic variations in our light curves. However, some dispersion remains due to 
photometric uncertainties, small differences in color and properties between the ensemble 
stars and the target objects, and the statistically correlated ("red" noise) that is present in 



all gro und-based transit surveys ( jPont et al.l l2006l : iBeatty fc Gaudil 120081 : Ivon Braun et al. 



2009b|) 



To estimate our remaining photometric uncertainties for every light curve, we calculated 
a weighted average magnitude and the standard deviation of the entire light curve about 
this average. Figure [7] shows the standard deviation of our BOKS light curves as a function 
of their r magnitude. As can be clearly seen, our brightest stars have a 4 milli-magnitude 
dispersion about the average magnitude, which sets our noise floor for this survey. For stars 
brighter than r ~ 14, saturation becomes an issue with seeing changes making the exact 
magnitude of saturation somewhat imprecise. 

As the apparent magnitude increases, the photometric errors also in crease. Plotted on 



Figur e [7| is an estimate of the photometric errors as a function of target flux ( lEverett fc Howell 
200 ih : 

where A^* is the number of ADUs from the star in the aperture, g is the gain of the CCD in 
electrons (2 e~ per ADU), Upi^ is the number of pixels in the aperture, Usky is the number 
of pixels in the annulus around the aperture used to measure the sky flux, Ngky is the flux 
in ADUs per pixel from the sky, and R is the rms read noise of the CCD in electrons (12 
e~^). For this comparison, we assumed the faintest value of the night sky brightness in our 
survey, which should give a lower bound to the true photometric uncertainties. As can be 
seen, this function is in good agreement with the lower edge of our error distribution. At the 
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faint end of our photometry, near r = 19, our uncertainties are ~50 milli-magnitudes per 
observation, substantially larger than the 10-30 milli-magnitude precisions needed to find 
transiting extrasolar planets. 

For completeness, we have also calculated the variability level expected from atmospheric 
scintillation alone (Young 1967; see also Young 1993a; Young 1993b; Badiah et al. 1994 for 
some important comments) in order to rule it out as a significant contributor of our highest 
precision photometric results. We find that the scintillation at unit airmass is approximately 
0.3 milli-magnitudes, about a factor of 15 lower than the photon noise from the brightest 
stars in our sample. This is an approximate value: Young (1993b) notes that the value 
can vary by up to a factor of two on a timescale of a few minutes. Nevertheless, given our 
measured photometric uncertainties, scintillation is not a major contributor to our error. 

5.2. Point Source Variability Statistics 

An important goal of BOKS is to determine the variability properties of stellar sources in 
general within one portion of the Kcplcr field. As mentioned previously, DAOFIND identified 
56,354 point sources in at least one BOKS survey image. Of these, 32,806 point sources were 
observed at least 1,000 times within our survey. Prom this subset, 25,861 point sources are 
located within 2.5 arcseconds of a source in the KIC78, and therefore we have additional 
photometric information. However, since the KIC78 data is non-contemporary with BOKS, 
any values derived from this comparison should be treated carefully. For statistical purposes, 
we excluded any source from the KIC78 that did not have a valid r magnitude, even if it 
had measured magnitudes in other bands. This left 22,340 sources for further analysis. 

We then determined a magnitude cut-off value for our statistical analysis of variability. 
Progressively fainter stars have larger photometric errors, are more likely to have poor quality 
light curves due to contamination from nearby bright stars, and in conditions of poor seeing, 
and may have incorrect rcccntcring by the aperture photometry algorithm, which can cause 
the aperture to recenter on a nearby brighter star. For these reasons, we restricted our 
variability statistics to the 13,786 stars brighter than r = 18.5, which also lie within 2.5"of 
a source within the KIC78 and which have over 1000 data points. We hereafter refer to this 
subsample as the BOKS-KIC sample. 

After applying the statistical tests from OPTICSTAT, we found 2,457 stars with > 5 
and r < 18.5 in the BOKS-KIC sample. We note that this number is approximate, as 
the reduced is strongly affected by a number of systematic uncertainties, such as residual 
variability in our ensemble stars, color mismatches between the ensemble and the target stars. 
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spatial structure in the extinction correction, and other effects. Additionally, photometric 
uncertainties and stellar variability applied to the KIC78 will move stars above and below 
the magnitude cutoff, creating an additional systematic uncertainty. 

The variability fraction found (~ 18%) can be compared to other variability surveys 
using identical search techniques. In an earlier study, Everett et al.l (120021 ) found a variability 
fraction of 17% over a five day survey period using a similar telescope and sampling rate, 
but using a less strict variabi l ity criteria (vf, > 3). In contrast, the varia bility study of 
NGC 2301 JXonrv et allboosl : iHowell et ahlboosl : ISukhbold fc Howell! boooh found a much 
larger variability fraction (56%) using an orthogonal transfer CCD observing mode over a 
12 night run, and with substantially better photometric precision 1.6 mmag). 



Tonry et al.l (120051 ) and iHowelll (120081 ) discussed a relation between the percentage of 
variable sources that will be found in any given photometric survey and that survey's pho- 
tometric uncertainty. From the NGC 2301 results, the cumulative fraction of variable stars 
found, as a function of quartile variability, x, (the quartile variability is ~ 1.5 times smaller 
than the standard deviation) is0 



f{<x) 



1.6 mmag 



X 



(2) 



For a best precision of 4 mmag, as we have in the BOKS survey, the percentage of variable 
sources would be expected to be ~40%, substantially larger than what was actually detected. 
Some of this difference may be due to d iffering stellar populations: NGC 2301 is a young 
operi cluster (250 Myr; Ixim et al.lboOlh . in contrast to NGC 6811 (575 Myr; ILuo et al. 
20091 ). and the level of stellar activity on the main sequence may be significantly different 
between the two clusters. The amount of background and foreground contamination may 
vary significantly between the two surveys as well. Given the position of the BOKS field, it 
is likely that many of the stars are field stars and therefore have lower amounts of variability. 

Of the 2,457 variable sources, 776 (32% of all the variables or f« 6% of the total BOKS- 
KIC sample) were found t o be periodic v ariabl e candidates (see §5.2.2p . This is significantly 
larger than the results o f Everett et al.l (120021 ). who found a 14% periodic fraction in their 



survey and iHowell et al.l (|2005[ ). who found a 13% periodic variable fraction in NGC 2301. 



We believe this is due to the substantially longer time coverage of BOKS, which should make 
it substantially easier for our search algorithm §5.2.21 to detect periodicity. The remaining 
variable stars appear to be non-periodic within the limits of our time sampling, photometric 
precision, and observational window. These non-periodic variable stars are presented in 



Again, we note that in this paper we use > 5 while the previous studies used > 3, which would 
allow more low amplitude variables into the samples. 
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Table 2, with the corresponding photometric information from the KIC given. We caution 
that the standard deviations given for these hght curves are hkely to be an overestimate: any 
photometric residuals in the ensemble stars and the effect of radiation events can increase 
this value significantly. 

There have been three other recent varia bility surveys of the regions of the BOKS field , 
all centered on the NGC 6811 open cluster (Ivan Cauteren et al.ll2005l : iRose fc Hintzl 120071 : 
Luo et al.l 120091 ). Unfortunately, the majority of the stars found to be variable by these 
su rveys are saturated on the BOKS images. The thre e exceptio n s, sta rs V8 and V9 found 
by Ivan Cauteren et al.l (120051 ) . and star V17 found by ILuo et al.l (120091 ) . were also detected 
by BOKS, and all three were found to be periodic variable stars. 



5.2.1. Non-variable Stars 

Our variability analysis found 10,881 stars brighter than r = 18.0 with no detected 
variability {xl <5). Figure [8] shows a color-magnitude diagram (CMD) of these stars. This 
CMD is typical of stellar fields within the galactic disk, showing a plume around a color 
of(7 — r = 0.4 — 0.6, corresponding to the main sequence turn-off for an old (> 10 Gyr) 
stellar population. The BOKS survey combined with the KIC does not go photometrically 
deep enough to detect the very faint and red low-mass stars in the disk population, which 
ty pically appear at q — r 1.4 and begin a.t r ^ 18.2. Similar CMDs are found in the work 



of Ide Jong et al.l ( l2010l . see their Figures 2 & 3) with data taken from the Sloan Extension 



for Galactic Understanding and Exploration (SEGUE). 

Figure ini shows a histogram of the number of photometrically constant stars as a function 
of r magnitude. The number of photometrically constant stars increases roughly linearly up 
to the final bin (r = 17.9). This reflects both the increase in the number of faint stars and 
our decreased sensitivity to variability for fainter stars. 

Some example light curves of stars that show little signs of variability over the timescale 
of our survey are given in Figure [TOl Note that in such long surveys, it is highly likely for 
a star to suffer at least one, and possibly more, hits from radiation events, such as cosmic 
rays. 



5.2.2. Periodic Variable Stars 



In or der to test each variable sta r light curve for per iodicity we applied the Lomb-Scargle 
method (ILombl Il976l : IScarglel Il982l ). as described by iPress et al.l (119921 ). The algorithm 
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produces a periodogram giving probabilities for sinusoidal signals in the light curves over 
a range of periods from our minimum sampled period of 1/60 days (0.4 hours), set by the 
spacing of consecutive exposures, up to 40 days, the maximum possible period spanned by 
the entire data set. We identified stars as possible periodic variables if: 1) OPTICSTAT 
returned a false probability of periodic variability of less than 1 x 10~^, 2) the amplitude 
was less than 2.5 magnitudes, and 3) the average magnitude was brighter than r=18.5. This 
yielded a large but manageable hst of candidates. We then visually inspected the light curves 
of these candidates after phasing them to the best-fit period. Only those light curves that 
had clear periodic signals, and whose light curves showed no sign of systematic effects were 
accepted as potential periodic variable stars. In order to reduce the effects of aliasing, we 
also removed from our sample any sources that had periods of 1 ± 0.025 days, though some 
of these objects may be genuine variable stars. 

We found that 776 stars from our variable sub-sample had periodic signals that ranged 
from ~0.2 day to ~33 days. The coordinates, mean r magnitudes, determined period and 
amplitude and the corresponding KIC information are presented in Table 3. Additionally, 
we have classified the periodic variables into approximate types, which are also presented 
in Table 3. Of these objects, 78 (~ 10%) show variability amplitudes larger than ~0.1 
magnitudes. Another 235 objects (~ 30%) have periodic amplitudes of 1-3% and periods of 
1-3 weeks, which is consistent with rotational modulation due to star spots. A significant 
number (A^ = 93; 12%) of the periodic variables remaining have periods less than 2 days 
and photometric amplitudes less than 0^.05. These shor t period low amplitude variables 
are likely to be pulsational variables such as 6 Scuti stars tereger 2000 ). 



We next compared the properties of the periodic variable stars as a function of period, 
stellar color, and magnitude, as these distributions give insight on stellar properties in gen- 
eral. It is well known that there are systematic changes in the fraction and am ount of stellar 



variability across the H-R diagram (lEyer fc Mowlavill2008l : ICiardi et al.ll201ll . and references 



therein), but the precise distributions are still under debate. Second, studies of these dis- 
tributions allow us to compare our results with other high precision variability surveys and 
provide confidence that the survey is valid. Figure [11] shows the overall distribution of peri- 
ods in our sample of periodic variables. The most notable features of this plot are the large 
number of stars with periods longer than 10 days and the peak in the period distribution 
near periods of ~1 day. Despite our attempts to remove them, this peak strongly suggests 
that a fraction of these variable stars have derived periods that are reflective of aliasing due 
to our observational window function, rather than their actual period. Additionally, the 
stars found to be periodic with the longest periods only have one or two measured periods, 
and the period-finding algorithm may have mistakenly flagged these objects, when in fact, 
they may not be periodic. More observations of these particular stars will be required to 
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fully address this issue. 

In contrast to the number of photometrically constant stars, a histogram of the number 
of periodic variables vs r magnitude, which is plotted in Figure [T21 shows a steep rise in 
the number of periodic variables between r ~ 14 and r ~ 16, which is followed by a steep 
decline in the number of periodic stars fainter than r ~ 17. This decline is due to the rapid 
loss of photometric sensitivity to low amplitude variations for the fainter stars. Figure [T3l 
shows a histogram of the number of periodic variables versus g — r color. The overall shape 
of this distribution is very similar to the color-magnitude distribution for constant stars (see 
Figure [8]). There may be a small excess of the bluest stars {g — r< 0.3) and the reddest stars 
{g — r> 1.0), but it is unclear whether this is statistically significant. In contrast. Figure [H] 
plots the period of periodic stars versus their g — r color. Most variables have colors that 
are similar to those of constant stars regardless of period. There seems to be a small excess 
of stars with < g — r < 0.5 among the periodic variables with the shortest periods (P^l 
day). 

Figure [T5] presents the distribution of amplitude for the periodic variables within the 



BOKS survey. As has been previously seen (IHowell et al.ll2005l ). the number of variable stars 



increases as the amplitude of variability d e creas es. The dashed line indicates the power law 



model of variability found by iTonry et al.l ( 120051 ). which has a slope of -2. The fit is in good 



agreement, giving further evidence to applicability of this model. 

Although a full accounting of the BOKS periodic variable stars would be soporific, we 
briefly discuss some of the more interesting stars here, and we discuss two extreme examples 
in § 15.41 Among the periodic stars with derived amplitudes larger than 10%, we find at 
least ten clear pulsational variables. An example light curve of this type of variable star is 
plotted in Figure dnj We have also detected two probable RR Lyrae stars within the field, 
with approximate periods of 0.53 and 0.56 days, respectively. A light curve of one of these 
objects is plotted in Figure [T71 

There are at least 32 eclipsing/contact binaries within the BOKS survey field, with 
periods varying from 0.13 to 6.10 days. Figures [TS] and ^\ give examples of two of these 
systems. From visual inspection of the light curves , the majority of the se systems are of the 



W Ursae Majoris subtype, as would be expected (IHoffman et al.l 120091 ). However, we have 
also found at least eight Algol-type binary systems. 
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5.3. Exoplanet Transit Candidates 

A primary goal of BOKS is to searcti for any signs of transiting extrasolar planets in 
the data set. To search the light curves for transits, we used a simple test to find and 
flag light curves containing at least one occurrence of a diminution in the relative flux with 
parameters specified below. The algorithm searched a light curve for time intervals when the 
mean magnitude is statistically significantly fainter than in the preceding and following data 
points, in other words, an inverse "top-hat" light curve. The algorithm stepped through 
time in each light curve testing a grid of possible interval starting and ending times and 
durations and reported back the most significant possible transit found for each light curve. 
A probability of significance is determined for a ll light curves for each possible transit by 



calculating a Student's t-test (jPress et al.l Il992l ) statistic comparing the mean magnitude 



during transit to the mean magnitude of combined pre- and post-transit data points. Those 
light curves with the most significant probabilities, specifically those that have a formal 
false probability less than 1 x 10"^^, were then subjected to further inspection. In order to 
avoid detecting too many spurious light curve fluctuations as transits, we next placed further 
restrictions on the events reported by the transit-finding algorithm. First, we searched for 
transit durations only between one and 12 hours. Second, we placed limits on the depth of 
transit that merited further study. Large planets orbiting all but the smallest M-dwarf stars 
will result in transit depths of 0.5 magnitudes or less. We therefore removed all possible 
transits that had depths larger than 0.5 magnitudes. 

For the purpose of completeness, we decided to use all of the stars in this analysis. It 
is extremely unlikely that faint stars will show a genuine transit event, but including these 
objects in the search allows us to test for other systematic effects in the algorithms. The 
entire BOKS sample contains 54,687 light curves that were all analyzed using OPTICSTAT. 
Using our detection limit described above, OPTICSTAT identified approximately 1,445 light 
curves with evidence of transit events. 

Each of these 1,445 light curves were then inspected by eye with careful attention given 
to additional criteria. We required at least two data points during transit and at least two 
data points in both the pre- and post-transit light curve, and we also required that the pre- 
and post-transit data points are separated by no more than 24 hours from the time of ingress 
or egress. The transit search algorithm used by OPTICSTAT searched for an inverted "top 
hat" shape in the light curve. However, there are a number of types of variability that can 
lead to a "top hat" shape that must be eliminated through a visual inspection of the light 
curve. Extrasolar planets will cause transits that have a fiat-bottomed appearance, so any 
sharp-bottomed transits were rejected from further consideration. It is not possible with our 
data set to observe secondary transits resulting from the planet passing behind its star, so 
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any light curve showing secondary transits was also removed from further study. Finally, 
in order to do follow-up observations we required that the light curve have at least two 
transit events. We require this to confirm the first transit-event and to obtain an accurate 
determination of the planet period before any follow-up observations are planned. 

The light curve shown in Figure [T9] was incorrectly identified by OPTICSTAT as a transit 
candidate, but was easily removed by our manual inspection criteria: the transits in this light 
curve are too deep (~ 0.6 magnitudes), sharp-bottomed, and there is an obvious secondary 
transit. Nearly all of the transit candidates detected by OPTICSTAT were rejected using 
the simple requirements we have outlined. 

At the end of our analysis we were left with three exoplanet candidates: BOKS-45069, 
BOKS-40959, and BOKS-52481. Some basic properties of these candidates are given in 
Table HI and a plot of their lightcurves is given in Figure [201 For each of the candidates, we 
then determined an approximate period by phasing the light curve to the best-fitting phase. 
In the case of BOKS-52481, we detected one full transit and only a portion of another 
transit, so the measured period is substantially less certain than the other two candidates. 
We also determined an approximate transit depth by averaging the closest 40 light curve 
points immediately before and after the transit to obtain a baseline. We present an example 
of each transit in Figure [211 

The properties derived are similar to other ground-based transit detections, but a de- 
tailed analysis, including follow-up spectroscopic and imaging observations of these candi- 
dates is presented in Howell et al. (2010). 

5.4. Two Specific Cases of Stellar Variability 

During any large area photometric survey, there is the potential of discovering unusual 
to rare objects. In the case of BOKS, we detail here two unusual variable stars that we have 
found in the survey. 

5.4.1. BOKS-459O6 

The first object, whose location is displayed in Figure[22l is known as BOKS-45906 (KIC 
9778689). On the first two clear nights of our survey (MJD 3980 & 3988), the star had a 
mean magnitude of r 20, though there were clear signs of variability of up to a magnitude 
in amplitude. The star continued to vary on both nightly and intra-nightly timescales. Then, 
between MJD 4004 and 4005, the star had an eruption, reaching a maximum of r = 16.6±0.01 
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on MJD 4006. It then declined in flux, returning to the approximate quiescent flux level on 
MJD 4009. The overall light curve is plotted in Figure [231 From the light curve, it is likely 
that this star is a newly discovered cataclysmic variable of the dwarf nova subtype. There 
i s scientific interest in the light curves of similar object s, both before and after eruption 
( lRobinsonlll975l : ICoUazzi et al.ll2009l : ISchaefer et al.ll2010[ ): therefore additional photometric 
monitoring of this source may be helpful. 



5.4.2. BOKS-53856 

The second interesting variable we have found is a blue star in the BOKS survey (a 
finding chart is displayed in Figure |2^ . known as BOKS-53856. From comparison to the 
KIC78, it has a measured color oi g — r = —0.46, making it the bluest stellar source in our 
field. Analysis of its light curve indicates periodic variability, with a period of 0.255 days, 
though of an unusual nature. The phased light curve is presented in Figure 123 

We obtained a 900 second spectrum of BOKS-53856 using the Kitt Peak 2.1-m tele- 
scope and the GoldCam spectrograph on UT 26 June 2008. We used the 300 1/mm grating 
(#32) with a one arcsecond slit to provide a mean spectral resolution of 2.4A across the 
full wavelength range. The spectra were reduced in the normal manner with observations 
of calibration lamps and spectrophotometric standard stars obtained before and after each 
sequence and bias and fiat frames collected in the afternoon. The final reduced spectrum is 
displayed in Figure [261 

The most obvious features in the spectrum are the blue continuum and strong Balmer 
lines, indicating a DA white dwarf type spectrum. In general, blue variables are either of 
low amplitude and consist of pulsations, or as in the case here, show larger variations and 
may be some sort of interacting binary. Holberg & Howell (2011; in preparation) present a 
further study of this star. 



6. Conclusions 

The goal of the BOKS survey was to constrain the amount and nature of variability in a 
subsection of NASA's Kepler mission field of view. The dedicated observations we conducted 
of the BOKS field allowed us to observe variability on various timescales from a few minutes 
to many days. The long-term observations also allowed us a reasonable opportunity to search 
for hot Jupiter exoplanet transits. 

Through our preliminary analysis of the variability in the BOKS field we have identified 
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~2,457 candidate variable stars with 776 candidate periodic variables. Most of this periodic 
variation can be attributed to rotating, low-mass stars with magnetic activity star/spots. 
We have also found over 90 6 Scuti stars, over 32 eclipsing binaries and contact binaries, as 
well as tens of large amplitude pulsators, such as RR Lyrae stars. Within the BOKS field of 
view, we have also identified at least three exoplanet candidates, all of which are undergoing 
observations by the Kepler mission. The comparison of ground-based and space-based transit 
observations should be beneficial to many future surveys. 

We would like to thank the entire staff of Case Western Reserve University Warner & 
Swassey observatory, including Heather L. Morrison, Charles Knox, and Colin Wallace for 
their invaluable assistance with the Burrell Schmidt. We are also extremely grateful to the 
observers of the AAVSO who performed many observations of our field concurrently. We 
thank Richard Wade for some useful discussions. We also thank the anonymous referee for 
several suggestions that improved the quality of this paper. This research has made use of 
the WEBDA database, operated at the Institute for Astronomy of the University of Vienna. 
This publication makes use of data products from the Two Micron All Sky Survey, which is 
a joint project of the University of Massachusetts and the Infrared Processing and Analysis 
Center/ California Institute of Technology, funded by the National Aeronautics and Space 
Administration and the National Science Foundation. 
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Table 1. Observing Log Summary 



Night 


JD 


Hours Used 


images 


Airmass Range 


Notes 


1 


2453980 


1.68 


25 


1.06 - 1.24 




2 


2453981 


0.00 






unusable due to weather 


3 


2453982 


0.00 






unusable due to weather 


4 


2453983 


A AA 

u.oo 






unusable due to weather 


5 


o A CT or\o A 

2453984 


A AA 

0.00 






unusable due to weather 


6 


2453985 


0.00 






unusable due to weather 


7 


2453986 


0.00 






1 1 1 J- -1-1 

unusable due to wcatner 


8 


O /I P" O A O T 

2453987 


0.00 






unusable due to weather 


9 


2453988 


0.42 


6 


1.19 - 1.25 




10 


C\ A P" A CI A 

2453989 


4.68 


56 


1.03 - 1.71 




11 


2453990 


0.00 






1 1 1 J- -1-1 

unusable due to weather 


12 


2453991 


0.00 






removed - CCD condensation 


13 


O /( P* O A AO 

2453992 


0.00 






removed - variable clouds 


14 


O /I CT O AAO 

2453993 


F A/^ 

5.96 


A 

69 


1.04 2.16 




15 


2453994 


F O T 

5.87 


O 1 

81 


1 n /I o no 

1.04 - 2.08 




16 


O A F O A A F 

2453995 


6.56 


95 


1.04 - 2.51 




17 


O /( F O AA/^ 

2453996 


6.33 


92 


1 AO A DO 

1.03 - 2.82 




18 


A /I P" O A A ^ 

2453997 


6.95 


98 


-i r\ A O A 1 

1.04 - 3.01 




19 


A /I F O AAO 

2453998 


A A A 

0.00 






removed - CCD electronics 


20 


O /( F O A A A 

2453999 


2.63 


38 


1.04 - 1.13 




21 


A /I F /I A A A 

2454000 


1.79 


24 


1.04 - 1.05 




22 


A /I F /I A A 1 

2454001 


6.54 


86 


1.04 - 2.55 




23 


A F A A A 

2454002 


0.00 






J "11 1 J 

removed - variable clouds 


24 


O /( F /( AAO 

2454003 


0.64 


A1 

91 


t n A O 

1.04 - 2.77 






/454004 


0.73 


no 
93 


1 n /I o Tn 
i.04 - ^.70 




26 


2454005 


6.60 


93 


1.04 - 2.72 




27 


2454006 


6.68 


96 


1.04 - 2.84 




28 


2454007 


5.86 


84 


1.04 - 2.08 




29 


2454008 


6.25 


90 


1.04 2.45 




30 


2454009 


4.75 


66 


1.04- 1.58 




31 


2454010 


5.76 


78 


1.03-2.18 




32 


2454011 


5.13 


57 


1.03 - 1.97 




33 


2454012 


5.28 


73 


1.03-2.21 




34 


2454013 


0.00 






unusable due to weather 


35 


2454014 


0.00 






unusable due to weather 


36 


2454015 


5.61 


74 


1.03 - 2.59 




37 


2454016 


0.00 






removed - lunar background 
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Table 1 — Continued 



Night JD Hours Used ^images Airmass Range Notes 

38 2454017 0.00 - unusable due to weather 

39 2454018 0.00 - unusable due to weather 

40 2454019 0.00 - unusable due to weather 



Table 4. BOKS Exoplanet Candidates 



BOKS 
ID 


KIC 
ID 


RA 
(J2000) 


Dec 
(J2000) 


r 

(mag) 


g _ j,a 

(mag) 


Period 
(days) 


Eclipse Depth 
(mag) 


Notes 


40959 


9595827 


19 39 27.667 


46 17 09.23 


15.1 


0.63 


3.9 


0.02 ± 0.01 




45069 


9838975 


19 40 08.003 


46 36 01.22 


16.1 


0.73 


2.6 


0.04 ± 0.01 




52481 


9597095 


19 41 18.802 


46 16 06.00 


15.9 


0.63 


7 


0.05 ± 0.01 


Period approximate 



Tor reference, the Sun is beheved to have a g — r color of 0.44 ± 0.02 (jBilir et al. 



2005 



Rodgers et al. 



20061) 
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Fig. 1. — The number of BOKS images taken as a function of survey time. Due to the 
effects of thick clouds, and occasional high winds, the numbers of nightly frames taken vary 
significantly throughout the survey. 
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Fig. 2. — An image of the BOKS field, created by combining tfie first 24 images of tfie 
survey. This image is approximately 101.' 5 by 49.' 5 in size. North is up and East is to the 
left. The open cluster NGC 6811, whose center is located at a = 19''37"*17^ S = +46'^23"*18^ 
(WEBDA database) is clearly visible in the center of the image. 
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Fig. 3. — The median seeing and median sky brightness for all of the r images in the BOKS 
survey. Filled squares denote the points used in the final light curve and open squares 
denote omitted points. Given the long length of the survey (40 nights), there are broad 
differences in these parameters. Note that the derived surface brightness comes from lunar 
phase differences, the presence and absence of clouds, and also the effect of astronomical 
twihght. The times of the full and new moon are shown for reference. The seeing values 
are measured from the images themselves, which have a pixel scale of 1.45", and are likely 
overestimates of the ambient seeing at the time of the observations. 
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Fig. 4. — The temporal window function for the BOKS data set. The sohd black line indicates 
the likelihood of viewing one transit event, the dashed red line indicates the likelihood of 
viewing two transit events, and the dotted blue line indicates the likelihood of viewing three 
events, as a function of the system's period. Given the ground-based nature of BOKS, there 
are strong aliasing effects at periods of integer days. Note that this window function is 
approximate: see §3 for discussion. 
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Fig. 5. — A plot of the magnitude correction derived from the ensemble stars for a single 
region over the entire BOKS survey period. Note that this correction takes into account 
the effects of airmass, seeing, and transparency simultaneously. The airmass effect of each 
night's observations is clearly visible in the data, as the BOKS field was close to zenith at 
the beginning of each night, and set as the night progressed. 
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Fig. 6. — The deviation of the BOKS astrometric coordinates to the coordinates found from 
the 2MASS survey. These deviations should be compared to the pixel scale of the Burrell, 
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Fig. 7. — The standard deviation of the hght curve of each star as a function of r magnitude, 
as determined by OPTICSTAT, for all stars brighter than an r magnitude of 20. Note that 
this plot includes every source, so variable stars will lie above the diagonal sequence that 
denotes the photometric error function. The upward bending of the points brighter than r 
= 14. is due to saturation effects. The thick diagonal hne denotes the expected photometric 
uncertainty for a source observed at new moon, and should denote the lower edge of the true 
uncertainty distribution. The solid horizontal line shows a standard deviation of 4mmag, the 
long horizontal line shows a standard deviation of lOmmag and the short dashed horizontal 
line shows a standard deviation of 50 mmag, for reference. 
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Fig. 8. — A g— r CMD for the 10,881 non-variable stars detected with this survey that have 
entries in the Kepler Input Catalog. Note the single broad sequence of stars at g — r ^ 0.5, 
which corresponds to stars at the main- sequence tur n-off for an old s tellar population. This 
CMD can be compared to the much deeper data of Ide Jong et al.l (120 lOl ). and shows that 
the BOKS field is typical of fields observed through the Milky Way's disk. 
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Fig. 9. — A comparison of the number of photometrically constant stars as a function of 
magnitude. As expected, the number of photometrically constant stars rises with apparent 
magnitude since the total number of stars increases, and our ability to search for photometric 
variations depends strongly on signal-to-noise. 
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Fig. 10. — An example of four stars with differing magnitudes that were found to be pho- 
tometrically constant throughout the survey, using the criteria for variability. From top 
to bottom the stars are BOKS-33934, BOKS-11505, BOKS-4700, and BOKS-3767. Note 
that the magnitude axis is substantially different for each star, ranging from 700 mmag for 
BOKS-33934 down to 50 mmag for BOKS-3767. Many of the most deviant points in each 
of the light curves are due to radiation events (cosmic rays) hitting the star's aperture. 
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Fig. 11. — A histogram of the number of periodic stars in our sample with r magnitudes 
brighter than 18.5 versus the logarithm of the derived period. Roughly half of the periodic 
stars we have identified have periods longer than 10 days. There is also a noticeable spike in 
the distribution near 1 day, suggesting that some of the periods we have identified are due 
to aliasing near the cadence of our observations, despite the fact we removed objects with 
periods of 1 ± 0.025 days (§5.2.2). 



-37- 



200 



150 



I 100 



50 



14 



15 



16 



17 



18 



19 



r (mag) 



Fig. 12. — A histogram of the number of periodic variable stars as a function of r magnitude. 
The number of variable stars detected rises with increasing magnitude down to r~16. The 
number of detected variables drops rapidly between r~17 and r~19. This reflects the loss 
of sensitivity to low amplitude variables, such as stars with rotational modulation, as the 
photometric uncertainties increase. 
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Fig. 13. — A histogram of the number of periodic variables with r brighter than 18.5 vs g— r 
color. This distribution has a single peak near g— r=0.6 and a tail that extends to g— r~1.6. 
This is similar to the color-magnitude diagram of non- variable stars (Figure [8]), except that 
there may be a small excess of variables near g— r=0.2. 
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Fig. 14. — A plot of the period of the periodic variable stars versus their g— r color. As 
can be clearly seen, most variables have colors that are similar to those of photometrically 
constant stars at all periods. For the shortest period stars, P<1 day, there seems to be a 
small excess of stars with 0<g— r<0.5. 
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Fig. 15. — The distribution of periodic variables found in BOKS, as a fun ction of phot o metri c 
amphtude. The dashed hne is the predicted function of variabihty from iTonry et al.l ( 120051 ). 
The error bars are the Poisson uncertainties for each bin. 
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Fig. 16. — This phased hght curve is an example of a pulsating variable star from our survey. 
The black points are the individual photometric measurements (and error bars). The blue 
squares are weighted averages calculated for every 26 points of the phased light curve. The 
red horizontal line at r ^ 15.18 is the weighted average r magnitude of the entire light curve. 
The green horizontal line at r 15.26 is the magnitude listed for the star in the Kepler 
Input Catalog (KIC). 
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Fig. 17. — This phased hght curve is one example of a RR Lyrae star found within the BOKS 
field. From visual inspection, the star appears to be of the RRab subtype. The symbols in 
this plot have identical properties to those found in Figure UHl except that only the mean 
magnitude is displayed horizontal line. 
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Fig. 18. — This phased hght curve is an example of an echpsing binary of the Algol type. 
The symbols are the same as in Figure [T71 
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Fig. 19. — This is an example of an Algol-type eclipsing variable star that passed the initial 
OPTICSTAT tests for a transiting extrasolar planet, but was removed after visual inspection. 
The very deep (0.6 magnitude) and V-shaped primary eclipse and the clear presence of a 
secondary eclipse rules this object out as a planetary transit system. The symbols are the 
same as in Figure [T71 
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Fig. 20. — The three highest quahty extrasolar candidates from the BOKS survey. In order 
from top to bottom, they are BOKS-45069, BOKS-40959, and BOKS-52481. The errorbars 
are omitted for clarity, but typical 2a uncertainties for these three light curves are 0.012, 
0.014, 0.012 magnitudes respectively. 
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Fig. 21. — The lightcurves of our three highest quaUty extrasolar candidates near the tran- 
siting phase. In order from top to bottom, they are BOKS-45069, BOKS-40959, and BOKS- 
52481. The errorbars are omitted for clarity, but typical 2a uncertainties for these three 
light curves are 0.012, 0.014, 0.012 magnitudes respectively. 
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Fig. 22. — The finding chart of a dwarf nova found by the BOKS survey. This finding chart 
is 5'x 5' in size, and North is up and East is to the left on this image. The object, known 
as BOKS 45906, is located at a = 19''40™16.22^ 6 = +46'^32'^48.23", and is centered in 
this image, surrounded by a circle for reference. The emission north by northwest of this 
candidate is a Schmidt ghost, and is unrelated to the source. See §5.41 for a discussion of this 
object. 




Fig. 23. — The light curve of a highly variable star found by the BOKS survey. This star, 
BOKS-45906, underwent a large photometric outburst, which lasted approximately five days. 
This behavior is typical of dwarf novae (see §5.41 for discussion). 
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Fig. 24. — The finding chart for a very blue object found in the BOKS survey. This finding 
chart is 5'x 5' in size, and North is up and East is to the left on this image. The star, known 
as BOKS-53856, is located at a = 19'*41"^31.35^ 6 = +46^^06^11.16" and is centered in this 
image, surrounded by a circle for reference. See §5.41 for a discussion of this object. 
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Fig. 25. — The phased hght curve of an extremely blue [g — r = -0.46) object in the BOKS 
field. The symbols are identical to Figure [T71 The light curve shows signs of periodic 
behavior, but with unusual structure. See §5.41 for further discussion. 
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Fig. 26. — A spectrum of BOKS-53856, a very blue object in the BOKS survey, obtained 
with the KPNO 2.1m telescope. Broad Balmer line features and a blue continuum are clearly 
visible, suggesting that this object is a white dwarf star. See §5.41 for further discussion. 
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